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ABSTRACT 

If the binding energy of the pulsar's surface is not so high (the case of a neutron 
star), both the negative and positive charges will flow out freely from the surface 
of the star. The annular free flow model for 7-ray emission of pulsars is suggested 
in this paper. It is emphasized that: (1). Two kinds of acceleration regions 
(annular and core) need to be taken into account. The annular acceleration 
region is defined by the magnetic field lines that cross the null charge surface 
within the light cylinder. (2). If the potential drop in the annular region of a 
pulsar is high enough (normally the cases of young pulsars), charges in both the 
annular and the core regions could be accelerated and produce primary gamma- 
rays. Secondary pairs are generated in both regions and stream outwards to 
power the broadband radiations. (3). The potential drop in the annular region 
grows more rapidly than that in the core region. The annular acceleration process 
is a key point to produce wide emission beams as observed. (4). The advantages 
of both the polar cap and outer gap models are retained in this model. The 
geometric properties of the 7-ray emission from the annular flow is analogous to 
that presented in a previous work by Qiao et al., which match the observations 
well. (5). Since charges with different signs leave the pulsar through the annular 
and the core regions, respectively, the current closure problem can be partially 
solved. 

Subject headings: pulsars: general — radiation mechanisms: non-thermal — 
stars: neutron — elementary particles 



1. Introduction 



After more than thirty years of intense study, the origin of the pulsed high energy emis- 
sion from rotation powered pulsars is still an unsolved problem. There are two types of radi- 
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ation models, i.e. the outer gap model ( Krause-Polstorff. Michel 1985a| jbl: lcheng et al. 1986 ) 
and the magnetic polar cap model. Within the polar cap models, two sub-types of model ex- 
ist. The vacuum gap model is based upon the assumpt i on of strong ion binding e nergy on the 



neutro n / strange star surface (IRuderman. Sutherland! Il975l: iGil. Sendyk I2OOOI: IZhang et al. 



1997br). The space-charge-limited flow (SCLF) model (IMichel 



1974 



Fawley et al. 



1977 



Hardinglll98ll ) assumes a low ion binding energy. However in the prev i ous p olar cap models, 
only the core cap has been considered, for a critical review see iMichell (119821 ). It is suggested 
in this paper that two kinds of accelerators, in both the annular region and the core region, 
must be taken into account. The a nnular regio n is def ined by the magnetic field lines that 
cross the null charge surface(NCS). iQiao et al.l (l2004al ) emphasized the importance of this 
annular region within the vacuum gap model. The model is found suitable to reproduce the 
7-ray pulse profiles as well as the radio pulse profiles. In this paper, we focus on the free flow 
case. Two major assumptions have been adopted in this paper, which are the key points of 
the model. They are: 1). For an oblique rotator, particles in the magnetosphere flow out 
though the light cylinder and the returning currents are located in the annular region; 2). 
The plasma in the magnetosphere is not fully charge-separated. 

In §2, we discuss our geometric setting and introduce the two polar flow regions. Pair 
production and 7-ray emission are discussed in §3, which focus on the properties of the 
annular acceleration region. The conclusions are presented in §4 with some discussions. 



2. Annular and core acceleration regions 

The open field line region of the pulsar magnetosphere is divided into two parts by the 
critical field lines. The part that contains the magnetic axis is the core region, while the 
other part is the annular region. The pulsar polar cap is also correspondingly divided into 
the core and the annular polar cap regions (Fig. [1]). For an aligned rotator, the radius from 
the magnetic pole to the edge of the core cap is = (2/3)^/^i?(i?f2/c)^''^, while the radius 
from the pole to the outer edge of the annular cap , i.e. the radius of the whole polar cap. 



is Tp = R{VlR/ cY^"^ (IRuderman. Sutherland! Il975l ). Here R and Vt are the radius and the 



angular velocity of the star, respectively, and c is the speed of light. 

For easy discussion in the following, throughout the text we will focus on anti-parallel 
rotators, i.e. ■ /i < 0. The case of ■ /i > could be derived by reversing the signs, and 
all the conclusions in this paper remain valid. 

For n ■ /i < 0, the Goldreich- Julian charge density (pgj) is positive in the region enclosed 
by the null charge surface (NCS). For the core- region field lines, positive charges flow out 
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Fig. 1. — A schematic diagram for the annular and the core regions for an inchned rotator. Q 
and fi are the rotational and magnetic axis, respectively. The null change surface (NCS) is the 
surface where the magnetic field is perpendicular to the rotation axis (i.e. ■ B = 0). Line 
'CFL' and 'LOF' are the critical field line and the "last" open field line respectively. Charged 
particles with opposite signs leave off the annular and the core cap regions, respectively. The 
upper boundaries of the pair formation front are denoted by the dashed lines. The "flash" 
symbols indicate the locations with the maximal possibility to initiate the pair-production 
cascade. 
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through the light cyhnder. The supply of positive charges from the surface compensate 
the deficit from the light cylinder. To maintain the invariance of the total charge of the 
pulsar, one needs a current that carries positive charges back to the star, or a current that 
carries negative charges away from the star. For the fully charge s eparated magnetosphere, 
only one sign of charge presents at a given location (IMichell Il979l ). In the annular region, 
the negative charged particles at lower altitude can not pass the region w ith positive GJ 



charge density at higher altitud e. An outer gap then forms beyond the NCS (lHollowaylll973 
Krause-Polstorff. Michej Il985al jbl: ISmith et al l boOlh . 



However, resent p ulsar magnetosphere simulations pr esented by so n ae aut hors (ISpitkovsky 



20061 : iTimokhinI |2006| ) give different results from that of ISmith et al.l (120011 ). Furthermore, 
there is no a prior justification for a fully charge-separated magnetosphere. Even if the mag- 
netosphere is initially charge-separated, the pair plasma generated from the outer gap will 
soon fill the annular flux tube, resulting in a quasi-neutral plasma. On the other hand, if 
one dismisses the conjecture of a fully charge-separated magnetosphere, another natural way 
to maintain the charge conservation of the pulsar is simply by extracting negative charges 
directly from the annular polar cap region (for the validity of the picture, see §3.3 for details). 
In this paper, we explore such a possibility. It is found that the negative charges stripped 
off from the stellar surface are naturally accelerated in the annular region. This leads to 
a particle acceleration model that keeps the geometrical advantages of both the polar cap 
and the outer gap models, which is found suitable to explain the 7-ray emission data (Dyks, 
Rudak 2003, Qiao et al. 2004). 

Since the positive and the negative charges are accelerated from the core and the annular 
regions, respectivel y, the pa r allel e lectri c fields (E h ) in th e two regions are opposite, as has 
been discussed by ISturrockl (Il97ll ) and iHollowayl (119751 ). As a result, vanishes at the 
boundary (i.e. critical field lines) between the annular and the core regions. also vanishes 
along the closed field lines. Thus it is equal potential along the closed field lines and the 
critical field lines. The potential on both the closed field lines and the critical field lines 
should be equal to the value at infinity. We also assume E\\ = at the star surface. When 
taking into account the 7-B process for pair production, there exist two accelerators, one at 
the annular region, and another at the core region. Furthermore, the pair formation front 
also moves to further distances near the magnetic pole. 

One important issue is whether the secondary pairs can screen the E\\ developed in each 
region. In the core region, since the primary charge density p has the same sign as the pqj, the 
produced secondary pairs tend to get polarized in the acceleration electric field and screen the 
field. This usually happens esp ecially if the primary 7-rays are produced through curvature 
radiation (IHarding et al.ll2002l ). For the annular flow, on the other hand, since the primary 



- 5 - 



charge density p has the opposite sign with respect to p^j, one has oc (p — pcj) < 0. As 
a result, pairs can not screen the acceleration electric field globally. Secondary pairs will be 
accelerated in the residual electric field. As a result, the charge acceleration region extends 
from the polar cap to higher locations in the magnetosphere. This potentially matches the 
geometric model proposed in Qiao et al. (2004). In the following, we will elaborate the idea 
more quantitatively. 



3. Acceleration, pair production and 7-ray emission 

3.1. Primary particle acceleration 

In fiat space-time, the calculation involving the 1-D Poisson's equati on and t he kin etic 
equation for charges to get the polar gap potential drop was performed by iMichell (119741 ). In 
the Kerr space-time with small a Lens e-Thirring angular velocity, the Poisson's equation is 
teeskinlEoOol : iMushmov. Tsvganlll992 L 



V- V-^*) = -47r(p-pGj), 



K 



(1) 



where k = 1 — rg/r, and Vg = 2GM/c^ is the gravitational radius. To the lowest order 
approximation, one has 

{VI - ujist) ■ B 

PGJ ^ TTT > (2) 

where ujuy = O.lbQR^r^^I^^ is the local Lense-Thirring angular velocity, and I45 is the 
moment of inertia of the star in unit of lO^^gcm^. 

In order to reveal the qualitative difference between the annular and the core cap regions 
(only within the primary accelerator), for simplicity we only solve the 1-dimensional Poisson's 
equation 

-47r(p-pGj). (3) 



K ds 



d ^ 
ds 

This equation is only valid when s < rp, but in Fig. [2] a tighter limit is placed. A full 3-D 
treatment is desirable to fully describe the electrodynamics of the system, and we postpone 
it to a future work. The 3-D calculation will be reduced to 1-D result, when the solution is 
confined in a region the transverse size of which is much larger than the longitudinal size, 
i.e. = 9V<9strans + 9'^ / ds\ong - d'^/dstrs.ns foY Strans >> Siong- Neglecting miuor energy 
loss due to radiation, the energy conservation law for particles in a given magnetic fiux tube 
is jmoc^ + = uiqc^, where we have used the condition $ = and E\\ = at the stellar 
surface. It should be noted that here the condition $ = at the star surface can only be used 
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in the 1-D calculation and should not be regarded as the boundary condition, since setting 
any $ at the star surface do not change the physical picture. The real physical boundary 
condition is E\\ = 0. This gives rriQC^d'y/ds = —qd^/ds. Here 7 is the Lorentz factor for 
the particles; s is the distance along the field line from the stellar surface; and m^^q are 
the mass and the charge of the particles we concern. In the annular region, the primary 
particles are electrons, so mo and q are the mass and the charge of the electrons. In the 
core region, mg and q are the mass and charge of the ions pulled out from the surface. In 
curved space time in the pulsar vicinity, the current conservation law can be expressed as 
^JT^J / B = ^JT^pv/B = const (for flat space time case it is pf/-B = const), where we assume 
that the charged particles do not cross magnetic field lines. Here v is the velocity for the 
particles, and J = pv is the current density. The charge density at a given height can be 
then expressed as p = ^/KQVoBpo/ y/KvBo, where the subscript '0' denotes the values at the 
stellar surface. Submitting the expression of p and p^j into eq. [31 we get 



d ^ 
ds 



1 d 



7) 



K ds 



X 



BvQ^fK^ Pgj 



BQV^/n pgjo. 



TRqC^ 



(4) 



where A is the reduced Debye wave length for the surface plasma, and a sign parameter 
X is introduced. For f2 ■ B < 0, one has x = ~1 for the core and the annular regions, 
respectively. For the ■ B > case, the following calculations are also applicable. The 
boundary condition at the stellar surface is p = p^ , in the core region, while p = — Pqj in the 
annular region. The later is based on the consideration that the negative polar cap current 
would eventually compensate the GJ current loss at the light cylinder. 

Equations ([3]) and (jl]) are solved numerically for both the core cap (CC) and the annu- 
lar cap (AC) regions, as presented in Fig. [2J Analytical approximate solutions are also 
derived. For the core cap, this is 7 ~ 1 + v^s/Acore, while for the annular cap it is 
7 ~ 1 + s^/A^^jj^ig^j.. We see that while the particle Lorentz factor increases linearly in 
the core re gion, it increases quadratical ly in the annular region, analogous to the vacuum 
gap model (IRuderman. Sutherlandlll975l ). The analytical solutions are also plotted in Fig. [21 
which show very consistent results with the numerical solutions. Notice that in order to 
compare the difference between the core and the annular regions, we have assumed that the 
accelerated particles in the core region are positrons. More realistic models involve positively 
charged ions. This would change the Lorentz factor calculation (lower panel of Fig. [2]) by a 
factor of mass ratio between the positron and the ion. 
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Fig. 2. — The accelerated potential drop ($) and Lorentz factor (7) as a function of distance 
s from stellar surface in the AC and CC for an aligned anti-parallel rotator, s is in units of cm 
and $ is in units of Volt. The legend "AC" and "CC" denote "annular cap" and "core cap", 
respectively. The subscripts 'Num' and 'Ana' indicate the results from numerical solution 
and analytic approximation, respectively. The shaded area is the schematic regions where 
the 1-D annular polar cap acceleration solution is no longer valid. The parameters used here 
are common values for 7-ray pulsar (i.e. p = 0.1 s, S = 10^^ G). 
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3.2. Location of the pair formation front 



The primary particles gain very high Lorentz factors (~ 10^ ~ 10^) within a short 
distance (~ 10^ ~ 10^ cm) (Fig. [2]). T hey will radiate 7-ray photons via curvature radi a- 



tion and inverse Compton scatte ring (IRuderman. Sutherland! 1 19751 : IZhang et al 



1997al lbl: 



Hardingl Il98ll : IZhang et al.l |2000| ). The energetic 7-ray photons will be converted into 
electron and positron pairs via the 7-B process. The condition to generate the secon- 
da ries is ■E'mrv-Bii2/44 > 1/ 15 if the magnetic fields are not close to the critical values 



g^7^^; for resonant inverse Compton scattering, the typi cal photon energy i s 
2h'jUB, where ujb = is the cyclotron frequency of electr on (IZhang et al.lll997al). 



(IRuderman. Sutherland! Il975l ). For curvature radiation, the typical 7-ray photon energy 

is —2c 
^■^ -'-^cur 

-E'rics = 

The naagnetic fields near t he stellar surface may not be pure dipolar (IRuderman. Sutherland 
I975I : iGil. Melikidzd I2OO2I ). We then generically assume that the curvature radius of the 
magnetic field lines is TZ, so that B_i = Bs/TZ, where B is the total magnetic field intensity. 
One can then estimate the typical length scale of pair production. The typical height of 
pair formation front (^o r gap h eight) in the c ore region h as been estimated prev iously, e.g. 
Ruderman. Sutherland! Jl975h : Izhang et all Jl997al . [200o[ ): Iflarding et aP J2002h . Here we 
focus on the annular region. In the annular cap, on e has 7 ^ 1 4 - ~ 0.21s^Bi2/P, 

where P is the pulsar period. Following the method of IZhang et al.l (120001 ). the height of the 
pair formation front in the annular region can be estimated as 



hmcs 



2.1 X lO^(cm) B^2"/'P'/'7^f 



(5) 



where the subscripts 'Cur' and 'RICS' denote the cases of curvature radiation, and resonant 
inverse Compton scattering, respectively. We can see that for typical parameters, the pair 
formation length is in the range of 10^ — 10'^ cm. Since the solutions for eqs.([3]) and (jl]) are 
valid only for acceleration of primary particles, the solution breaks above the pair formation 
front. In Fig. [2] we use a shaded region on the right to denote the region above the pair 
formation front. 



3.3. Secondary pairs 



The behavior of the secondary pairs in the core region has been widely discussed. The 
polarization of the pairs in the initial electric field generates an electric field with the opposite 
direction to compensate the initi al electric field. In m ost cases, the pair density is so high 
that the electric field is screened ( iHarding et al.ll2002l ). 



The case of the annular region is different. Whether the particles can be accelerated 



- 9 - 



depends on the electric field direction. If the electric field points towards the star in the 
annular region, the negative charges will be accelerated outwards and positive charges will 
be accelerated inwards. The only possibility for the electric field to vanish along the magnetic 
tube is p = pgj associated with trivial boundary condition. But the net free flow from the 
stellar surface is negative, which is opposite to the local p^j, thus the electric field can not 
be globally screened, although it may be partially screened locally and temporarily. The 
physical picture is that the negative charges move outward on a inward-pointing electric 
field background (with positive charges moving inward). The inward-pointing electric field 
is not a new idea. This is in fact discovered very long time ago by Holloway (1973), Cheng , 
Ruderman(1986), Michel(1979). They proposed the outer gap model, in which the positive 
charges move inward, which indicates the inward-pointing electric field or at least zero electric 
field. Since our basic assumption is a quasi-neutral plasma, the negative charges move 
outward, while positive charges move inward. 

An interesting possibility wit hin such a pictur e is that the annular flow is non-stationary, 
as has been discussed earlier by ISturrockl (Il97ll ). This is mainly caused by the interplay 



between the un-screened electric field and the screening fi e ld of the pairs. This may result 
in oscillating electric fields as envisaged by Leyinson et al.l ( 20051) , and may provide natural 



mechanisms to induce pulsar inward emission (jPyks et al. 




2005 



3.4. 7-ray luminosity 



The unscreened field extends at least to the NCS. Although with the existence of pairs 
E\\ is hard to model, electrons keep accelerating in this extended annular region, and almost 
achieve the maximum potential drop in the annular cap region, e.g. 



Max,Ann 



cP^ 



3/4^ 



+ 



a 



2 63/4 



(6) 



For a = 0, it is reduced to the value of the aligned rotator ( IRuderman. Sutherland! Il975l ). 
The maximum particle luminosity from the annular region is E'ann = ^PGjc$Max,Ann 
^2p-4(^2 + a), where A is the area of the annular cap region. We can see that the annular 
cap power is somewhat favorable to stars with large inclination angles. In Fig. [21 eq.([6]) is 
placed as the upper limit of $ (denoted by the shaded region above). 

An interesting inference from the above picture is that 7-rays are preferably emitted in 
the annular flow. Qiao et al. (2004) have investigated such a geometry model, and suggested 
that the model could well fit the observed wide-beam 7-ray pulse profiles as well as their 
radio pulse profiles. The acceleration model proposed in this paper provides the physical 
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basis of that geometric model. 

4. Conclusion and discussion 

The annular free flow is taken into account in this paper. In the calculation, we assume 
that E\\ = at the star surface and release the assumption that a pulsar magnetosphere 
must be charge separated. Several conclusions are derived. (1) There are two separate polar 
cap regions that accelerate particles with opposite signs; (2) Particle acceleration is more 
facilitated in the annular cap region thanks to its quadratic growth of the electric potential. 
It is then a favorable source of pair production. (3) Secondary pairs can not screen the 
parallel electric field in the annular region, so that electrons can be rc-accclcrated. The 
acceleration flow may be unsteady. (4) The annular acceleration region extends to the NCS 
or even beyond it. This leads to a fan-beam 7-ray emission, which has been found suitable to 
interpret the observed broad-band emission (Qiao et al. 2004). (5) Both the wide radiation 
beam observed in the high energy band and the current closure problem for pulsar can be 
addressed in the same time, if E\\ = at star surface is assumed. 

Whether the electric field is positive or negative in a given region not only depends on 
the sign of charges in the region but also the boundary conditions at the surface enclosing the 
region. Here the electric field is solved in a 1-D model. The potential at the star surface is 
arbitrarily chosen to be zero and the potential is not a boundary condition. This is not true 
for 3-D calculations, because the potential at the star surface and potential on the surface 
of tube is used as the boundary condition. In the SCLF case, boundary condition $ = is 
actually adopted at the star's surface, the charges can be accelerated in the charge-depleted 
region. However the boundary condition E\\ = is used at the star surface in our case, 
so charge excess rather than charge depletion is needed to accelerate particles in the core 
region. It should be noted that the Possion's equation is not well posed, if the electric field 
and the potential are both specified as the boundary conditions at the same places. 

The traditional outer gap model is based on the assumption of a fully charge-separated 
magnetosphere. Our model is based on the different assumption that negative charges are 
freely supplied from the surface of the pulsar. This leads to a different acceleration picture, 
in which the 7-ray emission region can be closer to the null surface. In principle, our model 
suggests that a vacuum outer gap should not exist. 

The re-acceleration and possible pair-production process beyond the pair formation front 
involves more complicated physics and is subject to further study. A more generalized 3-D 
treatment for both the core and annular accelerators is desirable. 
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Both the core and the annular accelerators produce pairs that can emit radio waves. 
The bi-drifting phenomenon (Champion et al. 2005) can be well interpreted by assuming 
that the observed radio emission comes from both regions (Qiao et al. 2004b). 



If pulsars are bare strange stars (IXu et al.lll999l ). vacuum gaps would be formed. The 



separation between the core and the annular region in that model has been discussed in Qiao 
et al. (2004). 

We are very grateful to the referee for valuable comments, and Dr. Han, J. L. and 
Mr. Zhu W.W. for useful discussions. This work is supported by NSF of China (10373002, 
10403001, 10273001). 
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